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D , ABSTRACT 
> '■ 

■ We give the detailed analyses for the gluonic-penguin effect on the Ktt and vrvr de- 

X 

■ cays of the B meson. In the standard model, it is shown that the ratio BR(B — > 
Ktt)/BR{B tttt) takes the value 0.5 ~ 3.0 with the strongly depending on the 
CP violating phase (p aiid the KM matrix element We obtain the constraint on 
the form factor by using the experimental branching ratio. It is also found that, in 
the two-Higgs-doublet model, the charged Higgs contribution which could enhance the 
B Xs'~f decay does not a sizable effect on our processes. The effect of the final state 
interaction on these processes is also discussed. 
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Until now, the rare B decays have been intensively studied in the standpoint of 
the standard model(SM) and also beyond the standard model. Especially, 6 — 57 
and b —>■ sg sub-processes have attracted one's attention in the circumstance that 
the experimental evidence has been found in CLEO These decays, induced by 
the flavour changing neutral current, are controlled by the one- loop penguin opera- 
tors which involve the important SM parameters such as the top-quark mass and the 
Kobayashi-Maskawa matrix elements Vts and Vf;,[0]. In our previous papers 0, we 
analyzed the inclusive decay B — > Xg'y and the exclusive decays B Kx'j, where 
Kx denotes the meson states in the sq{q = u or d) system, as well as the exclusive 
decays B — > Kx4> by including the nonstandard physical effects due to the charged 
Higgs contributions in the two-Higgs-doublet model(THDM) 

In this paper, we study the B Kit and B ^ mr decays, which are induced 
by both tree processes and gluonic penguin ones, in the SM model. Now we take 
into account QCD corrections, which were not included in the previous calculations [^], 
since the QCD corrections have been found to give the important contributions to the 
rare decays of the B meson. Recently, CLEO Collaboration reported the following 
experimental result 0: 

BR{B'^a + = (2.41!]:? ± 0-2) x 10"^ , (1) 

where the discrimination between the meson and the 7r+ meson has not been 
suceeded. We will predict the ratio of the B — > Kn decay width to the B ^ nn one, 
the ratio being almost independent of the form factor. The magnitude of the relevant 
form factor can be restricted by the experimental branching ratio of eq.(l) as shown 
later. Furtheremore, we will examine the nonstandard effects due to the charged Higgs 
contribution in the THDM. Finally, the effect of the phase shifts due to the final state 
interactions will be discussed as to our numerical results. 
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In the previous paper [Q, we give the formulation including the QCD corrections for 
the rare decays of the B meson. The four-quark operators and the magnetic-transition- 
type ones are given for the processes under consideration in the the Hamiltonian |^ , 

4Gf 



H. 



V2 



2 8 



i=l i=3 

where the factor t>g(q=u,t) is defined by 



(2) 



„, ^ J y,bV;, for 6 -> s 



V,bV;, forbad 
For the b s transition, each operator Oi is defined as follows 

O2 = {qLaYbLa){sLplt,qLI3) , 



03 = {SLal^bLa)i ^l' Lpll^lLfs) . 

5 quarks 

04 = (sL^rbLp){ E ^'lpIAJ , (4) 

5 quarks 
5 quarks 
5 quarks 

g 

Oj = -i^rnbSLa<y^''bRaqf,e„ , 
^8 = -^^rnbSLa(T^''T^pbRpq^et . 

For the b ^ d transition, the s-quark is replaced by the d-quark in each of eq.(4). 
In addition to these operators, we need a new operator Og, which is derived from 0% 
through the coupling of the virtual gluon to q'q': 

Og = imb^—qLo,(y^''T^pbRpq^qp-iuTpJ^ {q = s or d) , (5) 

where q^ denotes the four-momentum of the virtual gluons. The coefficients relevant 
to the processes B Kn and tttt are defined at the energy scale of mw asp|0] 

Ci{mw) = , C2{mw) = 1 , 



Csimw) = C^iniw) + 7 — ^ x 



x-f^ Xi In Xi . Xi , Xi — 6 ~t~ 2 



'2-^1 -xt il-xtf' 8^x^-1 (1-Xt 
C^5(mH.) = -^^G,{x,), (6) 

C4(miy) = C(,{mw) = ^'^^^ Gi{xt) , 

9o7r 

Crimw) = F{xt) , C'sl'Tily) = -gG'2(xt) . 

The functions Gi{xt) and F{xt) are given by 

- xj)(18 - Uxt - xj) - 2(4 - IQxt + Qx^^) Inxt 



Gi{xt) 

G2{Xt) = Xt 



'1 — Xt) (2 + 5xt — x^) + 6x4 In Xt 

ii-xtY 



nxt) = — -^[8x? + 5x,-7+^^;i^^^^lnx,], (7) 
24(1 -Xi)'^ {l-xt) 

with Xt = mf/m^r. 

In the following analysis, we take the value of the coefficient Cg(yu) being equal to 
GsifJ')- Although Cs{fi) does not include the full QCD correction of G'^{ji) in the leading 
log approximation, this replacement does not seriously affect the results numerically 
since the Og term is the next leading one compared to the operators O3, O4, O5 and Og. 
The similar operator O7, which is induced from O7 by the coupling of the virtual photon 
with g'g', is negligible due to a <^ as- We evolve the coefficients Ci{fi), by starting 
from the scale niw as given in eq.(6) to the scale fi = rrih = 4.58GeV, according to the 
renormalization group equation [|10|. Then, we obtain 



Ci{mb) = -0.240 , G2{mb) = 1.103 , 

Csirub) = 0.011 + 1.125C3(mH^) - 0.121Ci{mw) , 

Ci{mb) = -0.025 - 0.291C3(mvK) + 0.824:Ci{mw) , 

C5{mb) = 0.007 + 0.944C3(mvK) + 0.083C4(mvi/) , (8) 



Ceirrib) = -0.030 + 0.229C3{mw) + lAGbC^imw) , 

Cjirrib) = -0.199 + 0.629C3{mw) + 0.931C^{mw) + 0.675C7{mw) 

+ 0.091C8(mvK) , 

Csirrib) = -0.096 - 0.598C3{mw) + 1.029Ci{mw) + 0.709C8{mw) ■ 

However, the coefficients C^^mf,), 6*4(771;,), C^^mb) and C6(r7ib) do not completely involve 
the charm-quark loop effect. The charm-quark loop contributions are included by 



replacing these coefficients as follows |Tl 

C3(/i) - 
Cif,) - 

CM - 
CM - 

G{m^,q,IJ,) -- 



CM + '^^^G{m^,q, ^)C2{^l) , 
CM'^C{m,,q,fi)C2{fi) , 
CM + ^G{m,,q,l^)C2{l^) , 
Ceif^) - ^^4^C{mc,q,fi)C2{fi) , 

OTT 



4 / x(l — x) In 





777^ — — x) 



dx 



(9) 



where the parameter denotes the square of the four-momentum of the virtual gluons. 
These values of the coefficients are numerically given as follows: 



CMb) = -0.240, ^2(777,) 

Csirub) = 0.0152 - 0.00587 , CMb) = 

CMb) = 0.0118 - 0.0058Z , Ce{mb) = 

CMb) = -0.320, Csirub) 



= 1.103 , 

-0.0380 + 0.01747 
-0.0427 + 0.01747 
= -0.157 , 



(10) 



where is taken as ml/2 The imaginary part of these coefficients follows from 

the loop integral G{mc,q,f^). These values are in agreement with the ones given by 



Fleischer [11 



Let us begin with showing the decay amplitude of the B —>■ Kir process. By the 



use of the above operators, the decay amphtude is written as 

{K7r\H,ff\B) = 

+ 



AC 



V2 



(11) 



3,4,5,6,8' 

We use the factorization approximation in order to estimate the hadronic matrix ele- 
ment. This factorization assumption successfully works in the D meson and B meson 
decays [Q within the factor two as to the branching ratios. Under this assumption, 
the hadronic matrix elements of the above operators are given as 





7r^ 


\0i 


B',) = 


12^ 


|s7^75M|0)(7r- 1 117^6 15°) , 


{K^ 




\02 


B'a) = 


-i(A- 


S7^75w|0)(7r-|ll7^'6|5°), 


{K^ 


7r^ 




B'a) = 




|l7/.75n|0)(7r-|^l7^6|fi°). 


{K^ 






B'a) = 


-i(A- 


57^75M|0)(7r-|n7^'6|50), 


{K^ 




I05 


B',) = 


3^ 


Slur I 0)(7r- \uRbL \ 5°) , 


{K^ 






B',) = 


-2{K+ 1 


SlUr I 0) (vr^ 1 URbL B'^) , 


{K^ 


7r~ 


\0's 


B'a) = 


as 

rrih 

127r 


^q^[{K+ \s-f^-f,u\0){n- \ub 



+ {K+\sr,u\0){7r-\u^^b\B^,)) . 
By the use of these equations, the decay amplitude of -B — Kn is given as follows: 
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(12) 



+C4a^-(g^) + ^C5&^"(g2) + 2Ceb^^iQ'') + ^C8^c^"(g') ) ] , (13) 



a. 



37r 



where 

a^-(Q2) 
fe^'^(g') 
c^-(Q^) 



{K+\s^,^,u\Q){T,-\uYb\Bf) , 

(ir+|s75M|0)(7r-|n6|i?°), 

g'^ I S7^75M I 0)(7r- | ub \ B°) + {K' 



(14) 

S75n|0)(7r-|II7^6|B°)) . 



The hadronic matrix elements a^'^, b^'^ and c^^ are given in terms of the decay constant 
fx — 161MeV and the longitudinal form factor Fq'"{Q'^) as 



2\ 77B7r/'/^2^ 



m. 



{rris + mu){mb - my) 



m\ — m\ 



+ 



m 



K 



nil 



_2{mb-mu) nis + ruu 171% - ml 
where fx and the relevant form factors are defined by 

(X+ I S^n'JbU I 0) = fxQu. , 



(15) 



(16) 



(tt I ujf,b I 5^) = [pb+Ptt- 



m% — m 



with Q = Pb — P-K and then with = m\;. In the estimation of h^'^ and c^'", the 
equations of motion are used for quarks and anti-quarks. Then, the decay branching 
ratio is obtained by calculating 



BR{Bl K+TT-) = TB-^ I (X+TT- I //e// | B^,) 



?0\ 12 



'Sirm^ 



(17) 



The decay amplitude of i? — > tttt is given in the same way, 



(tt+tt- I H^ff I 5°) 



Gp 
V2 



+C4a'^'^(g^) + ^C,b''^Q'') + 2C,b''^Q') + ^Cs'^c''^Q'))] , (18) 



where 



a'^'^iQ') = /.(m^ - m^)Fo^'^(g^) , 



6""(g') 

c""(g') 



{md + mu){mb-mu) 



ml -ml 



+ 



mz 



m 



B 



2{mb - my) ma + mum%- ml 



a^-(g2) , 



with g2 ^ rnl and = 132MeV. 



(19) 
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In the calculation of eqs.(15)and (19), we used the following approximations: 

Ft (ml) ^ F,^-{ml) ^ F^iO) = Ff-{0) , (20) 
which are satisfied within the errors of a few percent in the pole dominance model of 



the form factor ||l3l. In our numerical calculations, we use the quark mass parameters 
as 



g rub = 4.58GeV or 5.12GeV, rric = 1.45GeV, = 160MeV, m„ = 5.7MeV and 
ma = 8.7MeV. 

Now we have left with only one unknown parameter Ft{0) in the B Kit 
and B ^ TiTT decay amplitudes except for the top-quark mass and the CP violating 
phase in the KM matrix. Then, we can predict the ratio of these decay branching 
ratios being independent of the form factor for the fixed phase 0, which is defined as 
Vub =1 Vub I exp(— i^). This ratio is almost free from the factorization assumption, 
because the ambiguity of this approximation cancels each other in the numerator and 
the denominator. We show the predicted ratio versus (p in the case of mi, = 4.58GeV 
and 5.12GeV with rrit = 140GeV in fig.l. Our result changes only 4% in the region 
mt = 120 ~ 180GeV. However, our result drastically depends on the value of | Vub/Vcb \ 
as shown in fig.l, where we use the experimental value | Vub/Vcb \= 0.08 ± 0.02[^ with 
\Vcb\= 0.045. 



fiff.l 



In order to know the contribution of the penguin process, we give the ratios of the 
penguin amplitudes to the tree amplitudes in the case of = 90°, nib = 4.58GeV and 
rrit = 140GeV for both B Kn and B ^ Tin processes as follows: 



A (penguin) 



y4(tree) 



4-22 x(^) for^^i^vr 
0.22 X i ) for B^nn 



(21] 



The penguin process dominates the B — ^ Ktt decay, but the tree process is not neg- 
hgible. On the other hand, the tree process gives main contributions to the B ^ Tin 
decay, ahhough the penguin process is still sizable. 

Since the experimentally observed branching ratio of i? — tttt + Ktt was given as 
shown in eq.(l), we can get the information of the form factor F^^'^^O) for the fixed 
I Vub/Vcb I and (f). We show the branching ratio versus F^'^{0) for = 90° and 30° in 
fig. 2. The experimental allowed region of the branching ratio is the one between the 
two horizontal dashed-lines in fig. 2. Then, we obtain F^'^[0) = 0.26 ~ 0.55, which is 
consistent with the one in the BSW model, F^^{0) = 0.33 



fig-2 

We comment on the contribution of the charged Higgs boson in the THDM as a 
typical new physics candidate. In contrast with the case of the B Xs'j decay, the 
charged Higgs contribution cannot provide so sizable enhancement on the B Ktt 
decays. This conclusion is in agreement with the one in the case oi B ^ Kx4' iii 
our previous paper ^ . Our predicted branching ratio increases only in the magnitude 
of around 10% for the case of itlh = 300GeVand cot f3 = 1, which follow from the 
experimental upper bound of the B Xsj decay |^. 

In the above analyses, we have neglected the final state interaction, which possibly 
affects the decay widths. The phase shifts due to the strong final state interactions 
have an effect on the magnitudes of the B Kit and B ^ im decay amplitudes. 
First decomposing their purely weak (this means "without the final state interaction") 
amplitudes according to the the final state iso-spins and introducing the corresponding 
strong phase factor for each iso-spin amplitude, we obtain the physical decay amplitudes 
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including the final state interaction. Then, we can readily rewrite, the physical decay 
amplitudes of 5 — > K^ti:~ in terms of the purely weak amplitudes of — > K^t:~ and 
5° ^ K^t:^ as 

^/2 

+ ^(e'^^- - l)(i^°7r° I H^ff I 5°)""="'= } , (22) 

where (^^tt = ^\/2 - (^3/2 and 

(tT+TT- I i/e// I S^)^'^^^ = e-^° {(tT+TT- I i/e// | 5^)"^"" 

+ ^(ei5.. _ ^)^^+^0 I ^^^^ I ^0^«,eafc | ^ ^23) 

where ^-^Tr = ~ ^2. However, the narrow resonances coupled to the Kn and tttt states 
are not expected at the energy scale of rri},. So, the large phase shifts are unlikely. 
Thus, our numerical results are not expected to be largely changed by the final state 
interaction. 

The summary is given as follows. We have studied the penguin effect of the B 
Kt: and 5 — > tttt decays considering the recent observed decay branching ratio by 
CLEO. We have predicted the ratio of these decay widths, which crucially depends on 
the CP violating phase (p and the still ambiguous KM matrix element | Vub \- The 
experimental information of the i^Tr/vrTr ratio will serve us these important parameters 
in the SM model. Furthermore, the determination of the form factor will test many 
models of B meson decays. We expect that the K-tt separation of the B meson decays 
will be done in the near future. 



10 



References 

[1] R. Ammar et al., CLEO Collaboration, Phys. Rev. Lett71(1993) 674. 

[2] M.Kobayashi and T.Maskawa, Prog. Theor. Phys. 49(1973) 652. 

[3] T.Hayashi, M.Matsuda and M.Tanimoto, Prog. Theor. Phys. 89(1993) 1047. 

[4] M.Tanimoto, K.Hirayama, T.Shinmoto and K.Senba, Z.Physik C48(1990) 99; 
A.J.Davies, G.C.Joshi and M.Matsuda, Phys. Rev. D44(1991) 2114: Z. Physik 
C52(1991) 97. 

[5] T.Hayashi, M.Matsuda and M.Tanimoto, preprint AUE-02-93(1993), submitted 
to Phys. Rev. D. 

[6] M.B.Gavela et al., Phys. Lett. 154B(1985) 425; 

L-L.Chau and H-Y.Cheng, Phys. Lett. 165B(1985) 429; 
M. Tanimoto, Phys. Lett. 218B(1989) 481. 

[7] M. Battle, et al., CLEO Collaboration, CLNS 93/1235, CLEO 93-13(1993). 

[8] B.Grinstein, R.Springer and M.B.Wise, Nucl. Phys. B268(1990) 269; 

R.Grigjanis, P. J.O'Donnell, M.Sutherland and H.Navelet, Phys. Lett. B286 
(1992) 413. 

[9] G.Buchalla, A.J.Buras and K.Harlander, Nucl. Phys. B337(1990) 313. 

[10] As to the renormalization group equation for the coefficients of OjS, see the refer- 
ences cited in Ref.[^. 

[11] R. Fleischer, Preprint of Technische Universitat Miinchen(TUM-T31-40/93). 



11 



[12] A.J. Buras, J.M. Gerard and R. Riickl, Nucl. Phys. B268(1986) 16. 
[13] M.Bauer, B. Stech and M. Wirbel, Z. Phys. C34(1987) 103. 
[14] Particle Data Group, Phys. Rev. D45(1992) II-l. 



12 



Figure Captions 



figure 1: The predicted ratios of BR{B^ K+7r-)/BR{B^ tt+tt") versus 
the phase for | Vub/Vcb |= 0.06,0.08 and 0.10. The sohd- and the dashed-hnes 
correspond to the predictions in the case of mi, = 4.58 and 5.12GeV with rrit = 140GeV, 
respectively. 

figure 2: The summed branching ratios of BR{B^ — > X+tt") and BR{B^ — > 
7r"'"7r~) versus the form factor Fq^'^(O) for | Vub/Vcb |= 0.06,0.08 and 0.10, being fixed 
777.6 — 4.58GeV and rrit — 140GeV. The sohd- and the dashed-hnes correspond to the 
predictions in the case oicf) — 90° and 30°, respectively. The upper bound and the lower 
one of the observed branching ratio are denoted by the two horizontal dashed-lines. 
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